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Afeatraot 

Two of the Bost challenging problems of Organo— 
aetallie chemistry (loosely defined) are pollntioa 
control with the large space velocities needed and nit- 
rogen fixation, a prooess so capably done by nature and 
so relatively poorly done by man (indnstry). For a 
computational chemist these problems are on the fringe 
of what is possible with conventional computers (large 
models needed and accurate energetics required). A 
summary of the algorithmic modification needed to 
address these problems on a vector processor such as the 
Cyber 205 and s sketch of our findings to date on deNOz 
catalysis and nitrogen fization are presented. 
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Introduction 


Two of the most chsllenging problems in Or g snom e t s 1 1 i c chem- 
istry (loosely defined) are pollution control with the large 
space velocities needed and nitrogen fixation, a process so 
capably done by nature and so relatively poorly done by nan 
(industry). For a computational chemist these problems (and 
other similar problems) are on the fringe of what is possible 
with conventional computers (large models needed and accurate 
energetics required). The advent of vector processors such as 
the Cyber 205 is making such studies feasible. A summary of the 
algorithmic modification needed to address these problems on a 
vector processor is presented in section I, a sketch of the 
findings to date for deNOz catalysis is presented in section 
II, and finally a sketch of the nitrogen fixation results is 
presented in section III. 

I. Algorithmic Modif ieat ioa. 

The advent of vector processors is leading to a reexamination 
of fundamental computational algorithms of general use to comp- 
utational chemists and the redesign of large scale codes. The 
present work illustrates both processes for the Cyber 205 comp- 
uter. Reexamination of fundamental algorithms is illustrated 
with an examination of the similarity transform, a matrix oper- 
ation of use to computational chemists. Large scale code rede- 
sign is examined through the implementation of a highly vec- 
torized HC-SCF code. 

A. Similarity Transform. A common sequence of matrix operations 
is the similarity transform 

C » A^ B A (1). 
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For compatat ional chemistr 7 applications the aatrices B and C are 
usually syaoetric and generally stored in lower diagonal forn. If 
the initial B matrix is expanded from upper diagonal form to full 
matrix representation vector operations are possible for both 
matrix multiplications. The linked triad instruction on the Cyber 
205 is utilized for the first matrix multiplication and a vector 
dot product operation is used for the second matrix multipli- 
cation. In principle one could transpose matrix A and to use the 
linked triad instruction for both matrix multiplications; 
however, in this case since we only want slightly more than half 
of the final results the vector dot product is preferable as it 
permits selective manipulation of the column indices I and J. As 
is apparent from Table I the vectorized matrix transformation 
represents a substantial improvement over scalar mode with 
enhancements ranging from a factor of 10 to a factor of 40. Note 
for the 300x300 matrix case we are still approximately a factor 
of 2 off the maximum rate for the Cyber 205. The consideration 
of an algorithm where several matrices are transformed at once is 
in order. In addition it should be noted from Table I that the 
expansion from lower diagonal form does not add a significant 
cost (less than 10 percent). Finally, it should be apparent that 
the MFLOPS rate will be independent of the number of orbitals 
involved (indices I and J); the vectorized loops run over number 
of functions not orbitals (indices K and L). 

B. SCF Coding Considerations. The fundamental kernel of self 
consistent field (SCF) codes in general ' is the energy 
expre s s ion 
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The integrals <X^lhlX^> and need only be 
evaluated once (for a given geometric point), stored conven- 
iently, and repeatively accessed during the orbital coefficient 
(C^) and density matrix element -(-Dj;D^i) optimization stages. 
For the Restricted Hartree Fock (RHF) wavefnnction dJ » 2, D^j = 
2, 1 * the remaining terms are zero.^ For wave- 
functions beyond RHF the wavefnnction optimization step repre- 
sents a vast majority of the time needed to variationally deter- 
mine E, that is, the calculation of the X^ integrals is usually 
relatively insignificant. For this reason initial ve c t or i zat ion 
efforts have concentrated on enhancing the time intensive stages 
of an MCSCF ( m u 1 1 i c on f i g ur a t i o n SCF) program. It is generally 
accepted that one of the most time intensive steps of a general 
MCSCF code is the 4 index transformation needed to convert the X^^ 
integrals to 0- integrals where 

0 . C ^X ( 5) . 
On scalar processors only the unique integrals are stored (the 
Canonical list) and the loops are structured so as to minimize 
the number of multiplications performed. On a vector processor 
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such as the Cyber 20S this step siaply amounts to two sequential 
applications of the matrix transformation described in (1). This 
transformation will proceed at vector speed provided that for a 
given ij pair all kl integrals are available fo'r k>l (this 
corresponds to an effective doubling of the integral file from 
its canonical length). This expansion of the canonical integral 
tape is accomplished through a straightforward two level bin sort 
written to take advantage of the 2 million 64 bit words available 
on the Cyber 205^. Since the vectorizable portions of this 
integral transform are contained in the matrix transform 
discussed above, the timing information in Table I applies here. 
Four index transformations for 50 basis functions will proceed at 
28 MFLOPS and 300 basis function transformations in general will 
achieve 82 MFLOPS. Enhancements over scalar computation on the 
Cyber 205 will range from a factor of 9 to a factor of 34 for 50 
to 300 basis functions. For example, a full integral 
transformation for 50 basis functions will maximumly take 28 
seconds and for 100 basis functions 10 minutes on the Cyber 205. 

For a wide class of useful w a ve f unc t i ons (open-shell HF and 
perfect p a i r i n g- g e n e r a 1 i z e d valence bond [GVB-PP] are two such 
examples) the one- and two— electron density matrices Dj and 
are expressible in diagonal form;^ that is, the only nonzero 
e lement s are 


Di 


2 f . , D ^ 4 

1 ' 1 j 


a-., and D^4 = b.. 

*ij* ij 
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The energy expresion (2) simplifies to 
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where 

= (ii/jj) and = (ij/ij) (8) 
are the nsnal Coulomb and exchange integrals. Restricting our 
attention to this class of wavefunction leads to particularly 
simple variational equations^ partitiouable into a step where 
occupied and virtual orbitals are mixed v a r i a t i o n a 1 1 y (OCBSE)^ 
and a step where independent occupied orbitals are mixed through 


pairwise rotations 


The OCBSE step utilizes terms 


representable as a vectorizable summation of and operators 


whe r e 




<X |J.|X^> = I Cj.ci (fivlon) 
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( 10 ) . 
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where a set of loops can be written (which are in linked triad 
form and will run at >170 MFLOPS for more than 50 basis 
functions) to evaluate the Ith hamiltonian (K runs from 1 to 


DO 300 J=1,NHAM 
A = A(I,J) 

B =■ Bd.J) 

DO 100 K = 1, MXS 
100 H(K)-H(K)+A*AJ(K. J) 
DO 200 K=1,MXS 
200 H(K)=>H(K)+B»AK(K, J) 


As the rotations step utilizes a subset of the above integrals. 


the needed v e c t o r i z a t i o n effort is narrowed down to rapidly 
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generating the teras in (10) 


If all 9,n teras a>n are stored 


for a given |iv the double sums in (10) can be reduced to a single 

dot product over a combined index y of length n(n-i’l)/2 

<X I J . lx > * 5 

|i ' 1 ' V L y y 

Y 

(13) 
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(14) 


Currently the are precalculated, stored, and used for an 
entire SCF iteration. Formulating the problem as in (13) permits 
vectors ranging from 1275 for 50 basis functions to 45150 for 300 
basis functions. This step will function at between 80 and 100 
MFLOPS representing enhancements of between 40 and 50 over scalar 
computation on the Cyber 205. Table II summarizes the timing for 
calculations ranging up to a 79 basis function calculation con- 
sisting of 4096 spatial configurations; that is, a GVB— PP( 1 2 /2 4 ) 
w a V e f un c t i o n . ^ If the calculation were stopped after the RHF 
step the SCF would represent less than 1% of the computational 
effort. Overall the GVB(12/24) wavefunction optimization repre- 
sents 14% of the total effort. This is in sharp contrast to 
computations on scalar computers where this step would account 
for greater than 95% of the effort. The timing for an SCF iter- 
ative cycle for three eases is broken down in Table III. Note 
that the time needed to generate the terms in (13) is comparable 
to that needed to diagonalize the variational hamiltonians 
(OCBSE) . 
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II. DeNOx Catalysis. 

Tlie catalytic redaction of nitrogen oxides has become increa- 
singly important in recent years due to legislation aimed at 
reducing emission levels from non-b i o log i ca 1 sources^. As Nitric 

7 

Oxide is the major NO^ component of exhaust streams research has 

focused on the reduction of nitric oxide. Both homogeneous and 

8 — 11 

heterogenous deNO^ studies have been performed . The use of 

base— metal catalysts is of particular interest due to their ready 

availability and low cost. A transition metal ion of singular 

importance in pollution control is Fe(II) either as the bulk 

oxide or ion exchanged into zeolites. These Iron systems have 

been demonstrated to catalyze the conversion of nitric oxide to 

8 9 

nitrogen with a co-reactant such as CO or H 2 ' . The mechanism 

1 2 

originally proposed by Shelef and Kumaer consists of a two 
stage oxidation reduction sequence. The initial step involves 
the coupling of two nitric oxides to fora nitrous oxide plus an 
Iron oxide. 


2N0 - ■■ ' ) NjO + '0' (15) 
The thus formed nitrous oxide is rapidly reduced by the cata- 
lyst®^ , 8d , 10 ^ 


NjO > Nj + ’0' (16) 

Completing the cycle the Iron oxide is reduced by reaction with 
carbon monoxide forming carbon dioxide plus the regenerated cata- 
lytic site. 

*0' + CO > CO 2 (17) 

Efforts have primarily been directed at characterizing reac- 
tion (15) as this is likely to be the kinetically most difficult 
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ft A 

step" . For homogeneous systems (15) has been suggested to 
involve an in t r a m o 1 a c n 1 a x coupling of nitrosyls to form a 
dinitrogen dioxide ligand^^* which rearranges to a bound cis 
hyponitr ite . 



1 2 3 

Metal hyponitrites have been established to either decompose to 
nitrous oxide and the metal oxide^^^ or react with carbon 
monoxide to from carbon dioxide and nitrous oxide^^^~°. 

It should be stressed that transition metal dinitrogen di- 
oxide complexes have never been isolated nor unambiguously 
detected. Further, only a single mononuclear transition metal 
hyponitrite complex has been identified 

In this section we report energetic support for the reaction 
sequence (18) for a model Fe(II) system: the dinitrosyl complex 
of Iron dichloride FeCl 2 (NO) 2 ^^. The relative energetics^^ and 
g e o m e t r i e s ^ for the chosen complex 1, its coupled cognate 
dinitrogen dioxide complex 2. and the cis hyponitrite product 3, 
are discussed below. We find that the coupled products are 
potentially accessible; 2 is only 29 kcal/mol higher in energy 
than 1 and 3 only another 19 kcal/mol higher. These species, 
though unobserved, should be viable given an appropriate ligand 
backbone. Addition of waters of hydration profoundly affects the 
relative energies of the hdyrated forms of 1. 2. and 3 (4. 5. and 
6 respectively). We find that Intermediates 5 and 6 are 
thermally accessible. Intermediate 5 is 24 kcal/mol more stable 
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than 4 and 6 is only 4 kcal/aol above 4 . This is not snprising 
as 1 is a Id electron systea. 2 is a 14 electron systea, and 3 is 
a 12 electron systea (unusual participation by the pi lone pairs 
vas not observed in the vavefunction of 3 or 6). 

A correlation of the bonding orbitals deaonstrates that the 

coupling reaction 1 to 2 or 4 to 5 will be thermally allowed 

(occupied reactant orbitals correlate with occupied product orb- 
1 7 

itals ). Further, the LUNO is a non-bonding d orbital of 
symaetry indicating that this correlation diagram will be valid 
for systems with up to 2 more electrons. Finally, one of the 
high lying occupied orbitals is a non— bonding d orbital 

suggesting that the correlation diagram will be valid for systems 
with up to two fewer electrons. Thus group VI through group VIII 
metal dications are potential active catalysts. 

Because Fe(II) dinitrosyls are structurally uncharacterized, 
because only a single transition aetal hyponitrite complex has 
been structurally characterized, and because dinitrogen dioxide 
complexes are unprecedented a detailed discussion of the bond 
distances and bond angles that were optimized is in order. We 
find the N-Fe-N angle for the dinitrosyl is 94.9 degrees, as 
expected for a {M(N0)2)® system^^^. The Fe-M distance of 1.69 A 
is in agreement with experimental structures for linear Iron 
dinitrosyls (1.66 to 1.71 A^®^). For the dinitrogen dioxide 

complex 2 we find a N-N distance of 1.53 A, longer than normal N- 
N single bonds (ranging from 1.402 A to 1.492 A^^) but still 
significantly shorter than that for free dinitrogen dioxide (2.24 
A ). This is consistent with substantial n i t r o g e n— n i t r o g e n sigma 
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bonding. Tie Fe~N distance found for the dinitrogen dioxide 
complex (2.23 A) is in accord with the Fe(II) nitrogen bond 
distance of 2.26 in [ F e ( C ^ H g NH) ^ ] [ F e ^ ( C 0 ) ^ 3 ] • Finally, for 

the eis hyponitrite complex 3 onr Fe-0 distance of 1.74 A 
compares favorably with 1.69 A (the sum of the ionic radii for 
0H~ (1.18 A) and an estimate for the ionic radius for four coor- 
dinate Fe(IV) (0.51 A)^^). Our N— N distance of 1.21 A is the same 
as the N-N distance determined by I— ray crystallography for 
[(Ph 2 P) 2 Pt(N 2 ® 2 ^^^^^* the only structurally chara c.t e r i z e d 
hypon i trite . 

Summarizing, we have demonstrated that (17) is a probable 
reaction sequence for group VI through group VIII transition 
metal deNO^ catalysts. Specifically our energetics and correla- 
tion diagram suggest that dinitrogen dioxides are thermo- 
dynamically and hine.tically accessible cognates of dinitrosyl 
complexes. Ve believe that these results can be extended to 
heterogeneous Fe(II) catalyzed deNO^ processes as well. In fact 
we speculate that the stretching frequencies observed by Hall 
at 1917 cm ^ and 1815 cm~^ are due to bound dinitrogen dioxide 

which is blue shifted relative to the free compound (which has 

2 3 — 1 — 1 

frequencies at 1870 cm and 1776 cm . Because the 

coordination sphere of Fe(II) ion exchanged into zeolites is 

thoughf* to contain three oxygen ligands our energetics suggest 

the frequencies assigned to a dinitrosyl are instead due to the 

hinetically accessible and thermodynamically favored dinitrogen 

dioxide moiety. Further, It should be noted that dinitrosyl 

stretching frequencies as high as 1900 ca~^ are rare. In 

conclusion we suggest that the kinetically (and thermodynamical- 
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ly) most difficult step is (17) is the isomerization of the 
dinitroges dioxide complex 2 (or 5) to the cis hyponitrite com- 
plex 3 (or 6 ). 

III. Nitrogen fixation. 

The fixation of dinitrogen is a reductive process of both bio- 
logical and large scale industrial interest. Thermodynamically 
the conversion of dinitrogen to ammonia is straightforward and 
the conversion to hydrazine is feasible under high pressures 
(A 62 pg for these processes are —7.9 kcal/mol and -*-22.0 kcal/mol 
respectively: if the pressure is increased to 100 atm then the 

^^298 hydrazine formation is -•■16.7 kcal/mol). 

In the known n i t r o g en— f i x ing organisms the catalytic reduction 

of dinitrogen is carried out by mo lybdoenzyme s known as nitro- 

genases . These nitrogen-fixing enzymes consist of two protein 

components, a Fe-No protein and a Pe protein. Further, an iron- 

molybdenum cofactor has been isolated from the Fe-Mo component 

protein of nitrogenase. In fact extracts of the Mo-Fe component 

from inactive mutant strains of microorganisms are activated by 

addition of this cofactor. Two models of the active site have 

been proposed that are consistent with Mossbauer and EPS spectro— 

scopic data- and EXAFS analysis of the Fo-Mo cofactor. 

Unfortunately the models of such active sites synthesized to date 

2 8 — 3 0 

do not reduce dinitrogen . 

Industrially, dinitrogen reduction occurs over an Iron cat- 
alyst at high temperatures and pressures. The rate determining 

3 1 

step is either the dissociative chem isorbt ion of dinitrogen 

2 * Nj > 2N-* (19) 
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or tlie siaple ehea i sorb t i on of an actiTated fora of dinitrogan 


• + Nj > Nj-* (20) 

Both of those processes are likely followed by rapid reaction 
with hydrogen (either aolecnlar hydrogen of eheaisorbed atomic 
hydrogen) . 

Thus, for both biological and industrial nitrification the 
activation of dinitrogen is a prerequisite for reaction with 
reductants such as hydrogen. Until very recently the observed 
forms of dinitrogen were bound to the metal with the nitrogen- 
nitrogen multiple bond largely intact ( no n- a c t i v a t e d ) . 

M.-=NjiN M=NaN.-H M j| (21) 

7 

Thus these model coapounds will only reduce dinitrogen under 

3 2 

rather harsh conditions . 

An understanding of a recently observed dinitrogen binding mode 
(analogous to organic azines) 

M«N-NaM (22) 

S 

will provide additional insight into biological and industrial 

nitrification. The reactivity and structural characteristics of 

3 3 

a new class of Tantalum complexes suggest the bonding pattern 

8 in (22). The Ta-N bond distances of 1.796 A and 1.840 A are 

quite similar to those observed in normal Tantalum imido 
3 3 

complexes (1.765 A to 1.77 A). In addition, reactions (23) and 

3 3 

(24) are both observed (reactions characteristic of metal- 
ligand multiple bonding). 
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M-0 


(23) 


M-N-R + R'jC-O > R'2C=*MR + 

M=N-N=M + 2 R 2 C-O )• R2C=»N-N=C2R + 2'M-O' (24) 

Finally, there is an observable 'activation' of the nitrogen- 
nitrogen bond (N-N bond distances of 1.282 A and 1.298 A compared 
to free dinitrogen which has a N— N bond distance of 1.0976 A). 

In this section we report energetic support for the kinetic 
and thermodynamic accessibility of 8 for molybdenum complexes. 
Our model consists of a bimetallic complex consisting of two 
Mo lybdenumte trachlor ide units bridged by a dinitrogen molecule. 
For this complex we have characterized the 'reaction path' 
connecting the two likely resonance structures 7 and 8 

Cl^Mo-K=t^.-MoCl^ 4 > Cl^Mo-N-N-MoCl^ (25) 

9 10 

We find local minima characteristic of each resonance structure 

indicating the 'resonance' interaction between these two forms is 

3 4 

not enough to result in a single averaged structure . However, 
the resonance interaction is sufficient to provide a very low 
barrier interconnecting them (less than 1 kcal/mol). Thermodyn- 
amically we find 9 to be 20 kcal/mol more stable than 10 for the 
tetrachloride ligand backbone. This thermodynamic difference 
could easily be overcome by an alteration of the ligand backbone 
and future studies will concentrate on this. Geometrically, for 
9 the Ho-N distance is 2.28 A and the N-N distance is 1.10 A and 
for 10 the Mo-N distance is 1.82 A and the N-N distance 1.23 A. 
This is in accord with a suggestion that the tetrachloride 
backbone does not fully activate the dinitrogen (a fully 
activated N— N distance should be on the order of 1.30 A). 
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Table I. Coaparisos of Scalar and Vector Matrix Transf oraat ions, 
(for Tarions sized astriees, tiaes in see.) 


Matrix 

Scalar 

(with 

Opt . ) 

Vector (tiaes x 100) 

Katio 

MFLOPS 

size 

First 

Se e ond 

Total 

Expand 

First 

Sec ond 

Total 

(S/V) 

( we 0 , ) 

NxN) 

Hnlt . 

Mult. 

Tiae 

Array 

Mnlt. 

Mnlt. 

Tiae 



50 

0.041 

0.083 

0.124 

0.053 

0.78 

0.51 

1.35 

9.1 

27.8 

100 

0.32 

0.55 

0.95 

0.23 

3.55 

2.59 

5.48 

14.8 

45.5 

150 

1.07 

2.58 

3.54 

0.51 

9.34 

5.91 

15.75 

21.7 

50.5 

200 

2.52 

5.74 

9.25 

1.01 

19.34 

14.32 

34.57 

25.7 

59.3 

250 

5.39 

14.35 

19.74 

1.83 

33 . 43 

25.54 

50.90 

32 . 4 

77 . 1 

300 

9.90 

27.14 

37.03 

2.92 

53.22 

42.23 

109.84 

33.7 

82.4 


Table II. Tiains Breabdovn for MC-SCF Energy Generation. 

(tiaes in seconds) 


Step 


Molecule/No. 

of basis functions 

H 2 O /7 FeClj 

•(H20)2/43FeCl 

2(N0)2/55 

FeCl2(NO)2(H20) 

Calculate 
One electron 
Integrals 

0 . 13 

35.4 

48 . 5 

81.0 

Calculate 
Two electron 
Integrals 

1.05 

85.5 

191.7 

535.5 

Sort Two 
Elec t r on 
Integrals 

0.05 

14.7 

94.3 

247.7 

Gene r a t e 
Extended Huckel 
Starting Guess 

_ ^ ^ 

0.8 

1.1 


Obtain 

Hartree Fock 
Energy 
(10 it . ) 

0 .11 

1.8 

3.1 

— 

Ob tain 
MC-SCF 
Energy 
(10 it.) 

— 

— 

72 . 5 

137 . 5 

Total Time 
% of Tiae 
HF 

MC-SCF 

1.35 

8.1 

140.3 

1.3 

411.2 

0 . 8 
17.5 

1001 . 7 
13 . 7 
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Table III. SCF Tiaing Breakdown for an Indiwidnal Cycle. 
(Tiaes in seconds, rates in MFLOPS) 


Wave f nnc t ion 
De script ion 

Generate 

and 

Matr ice s 

Transf orn 
J and 
Matrices 

OCBSE 

Orb it a 1 
Rotations 

Opt imize 
and b^j 

Total 


T ine 

Rate 

T iae 

Tiae 

Time 

Time 

Time 

H 2 O MBS HF 0 

.0001 

4.6 

0.006 

0.004 

— 

— 

0.011 

FeCl, • (H-O), 
HF 0 

.0082 

49.0 

0.017 

0.078 

— 

— 

0 .177 

FeCl- (N0)2 

HF 0 

.0310 

60.6 

0.034 

0.241 

— 

— 

0.306 

GVB(12/24) 2 

.012 

81.4 

2.832 

1.990 

0 .328 

0.091 

7.253 

FeCl2(NO)2(H2 
GVB(12/24) 4 

.302 

88.2 

5.322 

3 .515 

0.516 

0.090 

13.745 
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